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We present a general method to study weak-coupling instabilities of a large class of interacting electron 
models in a controlled and unbiased way. Quite generally, the electron gas is unstable towards a super- 
conducting state even in the absence of phonons, since high-energy spin fluctuations create an effective 
attraction between the quasi-particles. As an example, we show the occurrence of d-wave pairing in the 
repulsive Hubbard model in two dimensions. In one dimension or if the Fermi surface is nested, there are 
several competing instabilities. The required renormalization group formalism for this case is presented to 
lowest (one-loop) order on a most elementary level, connecting the idea of the "parquet summation" to the 
more modern concept of Wilson's effective action. The validity and restrictions of the one-loop approxi- 
mation are discussed in detail. As a result, three different renormalization group approaches known in the 
literature are shown to be equivalent within the regime of applicability. We also briefly discuss the open 
problem of a two-dimensional Fermi system at Van Hove filling without nesting. 
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1 Introduction 

The theoretical understanding of systems with strongly correlated electrons is one of the most important 
objectives of current research in physics and the key for elucidating a number of phenomena which have 
been observed in newly designed materials. The most spectacular example is obviously high-temperature 
superconductivity discovered 1986 in the cuprates, but intriguing questions have also been raised by other 
transition metal oxides, by organic compounds etc. 

Although the (bare) Coulomb interaction between electrons is generally strong in materials, it is im- 
portant to control at least the limit of weak interactions, which is already a difficult task. Furthermore, 
arbitrarily weak interactions lead to strong correlations at low temperature. It is thus worthwhile to study 
the weak-coupling problem thoroughly. 

As will be shown in detail, naive perturbation theory breaks down at low temperature and has to be 
improved by the so-called renormalization group (RG). RG concepts have been used in very different 
fields of modern physics and can have quite different meanings for different people. The notion was first 
introduced by Stiickelberg and Petermann 1 1 1 and independently by Gell-Mann and Low 1 2 1 in the context 
of quantum field theories (such as quantum electrodynamics) in order to cope with infinities that appear in 
naive perturbation theory. 

In the early 1970s, Kadanoff and Wilson have associated the RG to the procedure of mode elimination 
in classical statistical mechanics (or systems of bosons) OB). In Wilson's formulation, the RG trans- 
formation consists in integrating out some degrees of freedom of the system and including them in the 
renormalization of some parameters (for example coupling constants). This alternative formulation of the 
RG idea proved tremendously successful for the analysis of critical behavior in the vicinity of second order 
phase transitions Q. 

In general, a RG transformation is some change of the length or energy scale and the RG equations 
describe the response of the system as the length or energy scale is changed. In the interacting electron 
problem considered here, the energy scale is given by the temperature or alternatively by a cutoff A in the 
band energy of the electrons. The RG equations describe the change of a two-particle Green's function as 
a function of the energy scale A. In many cases, the RG flow to low energies produces a singularity at a 
finite energy scale A c . This is interpreted as a transition into a strongly correlated state. The energy scale 
A c is then comparable to the mean-field transition temperature towards an ordered state. 

It turns out that the one-loop (i.e. the leading order) approximation of the RG is equivalent to the 
so-called parquet approximation. This method was developed by the Soviet school 1 6 1 in order to treat 
different diverging fluctuations on an equal footing. It was successfully applied to one-dimensional con- 
ductors 1 7 8 9 1, to the Kondo problem 1 10 1 and to the A -ray absorption edge singularity problem in metals 
II 1111 21 . A detailed description of the method can be found in II II . In our derivation of the one-loop RG 
equations, we will follow mainly the parquet philosophy. 

The application of renormalization group (RG) ideas to fermionic problems in more than one dimension 
started in the mathematical physicists community 1131 [14 1 and has led to a considerable progress during 
the last decade 1 . For example, a series of rigorous studies has shown that the Landau Fermi liquid theory is 
stable for the D = 2 jellium model at not too low temperatures such that \U log T| < const, where U is the 
strength of the local interaction 1221 1231 . For lower temperatures, the properties of the interacting system 
are no longer analytically connected to those of the non-interacting system. The usual interpretation is a 
phase transition into a superconducting state, although there exists so far no rigorous statement about the 
regime T <T C . 



Pedagogical introductions can be found in [ 15 16 17]. A precise relationship between Fermi liquid theory and the renormal- 
ization group has been established in 1 1 8 19 20 2.1 J. 
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A number of different numerical schemes to compute the one-loop RG flow of D = 2 Fermion mod- 
els have been developed during the last years. They have revealed many interesting and appealing results 
such as the appearance of rf-wave superconductivity in the repulsive Hubbard model close to half filling 
1241 [25 1 , a spontaneous breaking of the lattice symmetry (Pomeranchuk instability) 1 26 1 and a phase with 
suppressed uniform spin and charge susceptibilities ("insulating spin liquid") 1 27 1 . The method was further 
developed to a finite temperature RG scheme in [28 1. An interesting extension of the RG formalism to ob- 
tain dynamical and non-equilibrium properties of certain impurity models has been presented in Ref. |29|. 

In this article, we have chosen a comparatively pedestrian approach to the weakly interacting electron 
problem and its RG treatment with the aim of being intuitive and didactical. The paper is organized as 
follows. In Section|2] we show how naive perturbation theory breaks down as a consequence of infrared 
divergences. In the case of a generic Fermi surface in more than one dimension, these divergences can be 
accounted for in a controlled way using the ladder approximation. This method is presented in Section|5J 
where we also discuss its relation to Kohn-Luttinger superconductivity and BCS-theory. This part is fol- 
lowed by a specific application in Section|4] showing the occurrence of superconductivity in the repulsive 
Hubbard model in two dimensions. 

The remaining part of this article addresses the more general situation, where there is more than one 
competing weak coupling instability. This is generally the case in one dimension and in more than one 
dimension for special (nested) Fermi surfaces, which are defined in Section [5] If the Fermi surface is 
nested, the ladder approximation is no longer justified and requires a non trivial extension. The necessary 
formalism is developed in Section [6] and We first derive several exact Bethe-Salpeter equations which 
are then used to establish the one-loop RG equations. The region of applicability of the one-loop RG 
is discussed in Section [8] In Section |9j we relate our approach to the Wilsonian RG concept and we 
finally compare several recent RG approaches known in the literature in Section^] In the Conclusion, we 
comment upon two-dimensional (2D) electron systems at the Van Hove filling. 

Some more technical aspects are treated in four appendices. Appendix lAl contains basic definitions 
and properties of the two-particle Green's function and the vertex function. In Appendix 151 we calculate 
a specific example of a third order parquet diagram, in order to demonstrate the general properties to 
which the main text refers. In Appendix Icl the RG formalism is extended to calculate linear response 
functions. Finally, some basic definitions and identities of the functional integral formulation are given in 
AppendixlDl 

2 Infrared divergences in naive perturbation theory 

We consider a general single-band model of interacting electrons on a D-dimensional periodic lattice. The 
Hamiltonian is of the form H = Hq + Hj, where 



is the most general two-body interaction respecting the global spin rotation symmetry, translational invari- 
ance and particle number conservation. V is the number of lattice sites (i.e. the system volume), cl and 
c ka are tne usua l creation and annihilation operators of electrons with lattice momentum k and spin index 
a =|, J,. The single-particle dispersion e and the coupling function g are smooth functions of the momenta, 
if the hopping matrix and the interactions are short-ranged in real space. 

Our aim is to calculate the imaginary time Green's functions. Special attention is given to the vertex 
function T(k%, ^3), which is defined as the one-particle irreducible part of the two-particle Green's 




(1) 



is the non-interacting part and 




(2) 
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function. The D+ 1 -dimensional frequency-momentum vectors ki = (fcoi,kj), involve the Matsubara 
frequencies fcoi which run over the odd multiples of tr/f3, where (3 is the inverse temperature. r(fci, fc 2 , fc 3 ) 
represents a scattering amplitude of two particles with incoming frequency-momentum and spin labels 
(fci, j) and (fc 2 , 1) into outgoing labels (fc 3 , J.) and (fci + fc 2 — fc 3 , j). A precise definition of the vertex 
function as well as some basic symmetry properties are given in AppendixlAl 

The vertex is always a regular function at high enough temperatures. At low temperatures however, 
the vertex can have poles. These poles are usually interpreted as the manifestation of an instability of the 
metallic state such as superconductivity or different kinds of density waves, since a divergence of the vertex 
implies a divergence of some generalized susceptibility. 

Since we are concerned with weak interactions, the first approach is standard perturbation theory in Hj. 
A calculation of the vertex function to second order in perturbation theory gives 

r(fci, fc 2 , fc 3 ) = -fftki.ka.ka) + PP + PHI + PH2 (3) 

PP = W 2 5(kl ' k2 ' k _ P) C ^ C ( fc - P)»(P> k ~ P> k 3)' 
p 

PH1 = ^F^ff(ki,P + qi,k 3 )C(p)C(p + (Zi) 5 (p,k 2 ,p + q 1 ), 



p 

PH2 = -L C (P) C (P + 92) [-2<?(k x , p, p + q 2 ) 5 (p + q 2 , k 2 , k 3 ) 

p 

+g(p, ki, p + q 2 )#(p + q 2 , k 2 , k 3 ) + g(ki, p, p + q 2 )#(p + Q2, k 2 , p)] , 

where k = k% + fc 2 , q% = fc 3 — fci, g 2 = fc 3 — fc 2 and the single-particle Green's function of the non- 
interacting system is given by 

C(k) = "77 ~ ~Z i (4) 

iko - 

where £k = — /i. The representation of the second order contributions in terms of Feynman diagrams is 
shown in Fig.^ The internal electron lines stand for free electron propagators C and the wavy interaction 
lines stand for coupling functions g, that depend on the incoming and outgoing momenta. The convention 
is that the spin index is conserved along the fermion lines. The minus sign in the first of the three PH2 
diagrams comes from the fermion loop and the factor 2 in the same diagram from the sum over the spin 
index in the fermion loop. PP is referred to as the particle-particle (p-p) diagram and PHI and PH2 as 
particle-hole (p-h) diagrams. 

Suppose for a moment that g is constant (= U in the Hubbard model). The p-p diagram is then given 
by g 2 multiplied by the p-p bubble 

p 

J_ "(gp) + "(£k- P ) - 1 (6) 

where n(£) = (1 + exp is the Fermi distribution function. This quantity is diverging at k — as the 
temperature goes to zero. At a finite temperature T = one finds 

f W , , tanh BL , s W 
BPP(0) = / d£ v(0 ——2- = KO) log — + "finite", (7) 



w 



2£ v ' T 
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Fig. 1 Second order diagrams for the one-particle irreducible vertex T(ki, ki, k^). k± is determined by energy- 
momentum conservation. 

where W is the bandwidth, v{u>) = 1/V J2k <^(£k — w) is the density of states, and "finite" is any contri- 
bution with a finite limit T — > 0. The only two assumptions are reflection symmetry £k = £_k and that the 
density of states z/(£) is an analytic function of £ at the Fermi level. Under the same assumption it can be 
shown that the largest term of order n in the coupling diverges as log™ -1 (W/T). 

It is clear that naive perturbation theory breaks down at large enough values of the coupling or low 
enough temperature, such that |<7log(M /r /T)| ~ 1, Alternatively, one can organize the perturbation series 
in the form 

r = ai g + a 2 g 2 H , (8) 

where 01,02,... are analytic functions of the product g \og(W/T). In the following we concentrate on the 
leading order contributions, i.e. we calculate only the first term in this expansion. This means that diagrams 
of a given order in g are calculated to leading logarithmic precision and all the subleading contributions 
such as the term "finite" in (7) are neglected. In other words, g is considered to be small, but g \og(W/T) 
is not. 

One can regularize the zero temperature perturbation theory by introducing artificially an infrared cutoff 
A. There are several ways to do this. One method of avoiding the singularity is to remove from the theory 
the single-particle quantum states which are too close to the Fermi surface. This amounts to replace the 
bare propagator by 

C A (k) = 6(|a| - A)C(fc), (9) 

where 8 is the Heavyside step function. In this case, the vertex function becomes cutoff-dependent and we 
use the notation Ta(A;i, &2, fe)- In the presence of the cutoff, the logarithms in Eqs. and l|8} are replaced 
by log ma) ,j^ A | ■ Alternatively, one could also use the quantity y^o~+~^k mst ead of to introduce the 
cutoff or replace O by a smooth function without changing the result to leading logarithmic order. We are 
thus free to calculate at zero temperature, introduce the cutoff in the most convenient way for calculations 
and in the end of the calculation replace A by T. It is clear that the correspondence A <-» T is only correct 
to leading logarithmic order. In a more elaborate calculation including subleading contributions, the result 
depends on how the cutoff is introduced. 
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1 




2 



3 



Fig. 2 The structure of diagrams, which are reducible in the p-p channel. The grey circles represent any sub-diagram. 



Fig. 3 The structure of p-p ladder diagrams. The wavy lines stand for — 1^ and the electron lines represent full 
propagators Ga, in the presence of the infrared cutoff A. 

3 Ladder approximation and BCS theory 

We now calculate the vertex to leading order in the expansion (|8). In a generic situation in more than one 
dimension without any fine tuned parameter, the only diverging terms occur in the p-p diagram. The other 
second order terms have a finite limit for T, A — * except in one dimension or for special situations, which 
will be discussed later. Having identified the p-p diagram as the main source of the infrared singularity, we 
define the class of p-p irreducible diagrams. A diagram contributing to the two-particle scattering vertex 
is called reducible in the p-p channel if it is of the form shown in Fig. [2] i.e. if it can be divided into two 
disconnected pieces by cutting two particle lines. In the opposite case it is called irreducible in the p-p 
channel or simply p-p irreducible. 

In the following, we use the notation I pp (ki : k%, k^) for the sum over all p-p irreducible diagrams. 
To second order in perturbation theory, this quantity is given by Eq. with the singular PP contribution 
omitted. It can thus be calculated by naive perturbation theory without the problem of infrared divergences. 

The exact full vertex is given in terms of the p-p irreducible part and the single-electron propagator by 
a summation over the so-called ladder diagrams (Fig. [3}- The analytical expression of this series is 




k 2 



P1-P2 



(10) 




In principle, the internal electron propagators in Fig.[3]are exact single-particle Green's functions G\(p) = 
0(|£p| — A)/ (ipo — £ p — Sa(p)), where Sa(p) is the self-energy in the presence of the cutoff. Self-energy 
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corrections have three main effects. First they change the shape and location of the Fermi surface, second 
they modify the properties of the single particle dispersion (the Fermi velocity) and third they lead to 
a reduction of the quasi-particle weight. The deformation of the Fermi surface requires in general the 
introduction of counter-terms 2 . We will ignore this effect and assume that in the weak coupling limit, 
the interacting Fermi surface is not so different from the non interacting one such that the difference is 
not crucial. Perturbative corrections to the Fermi velocity V p E(p) and the quasi-particle weight z = 
( 1 + id Po £ (p) ) ~ 1 are finite as T — > and therefore do not contribute to leading logarithmic order. We will 
therefore neglect self-energy corrections in the remaining part of this section. 

As argued above, the p-p irreducible vertex function can be calculated within naive perturbation 

theory. Let us thus focus on how to sum the infinite series Eq. (II Ot . for given J^ s . 

For simplicity, we set the frequency qo equal to zero and neglect the frequency-dependence of I BCS (fe, k'). 
This point will be justified later. We focus on small values of q and A, because this is where the infrared 
divergence appears. The p-p irreducible vertex has no infrared singularity and can thus safely be replaced 
by 

I BCS (k,h>)= lim I^ S (k,k'). (12) 

In this way, I BCS becomes a real-valued and symmetric function of the momenta k and k'. One can then 
define a Hermitian operator 3 , which acts on a function /(k) as 

/(k) — * FsfM ^ /BCS(k ' k,) ^ (k ' )/(k,) ' (13) 

where _D^ p q (k) = 1//3 53fc ^AqW ' s tne non-interacting frequency-summed pair propagator and i?^ p (q) 
is the p-p bubble in the presence of the cutoff. The factor l/B 1 ^ J (q) has been introduced in such a way, 
that Eq. (1131 acquires a finite limit for A, q — » 0. 

This implies that I BCS has a spectral representation given by 



I BCS (k, k') = ^ A„ e n (k)e„(k'), (14) 



where A„ are the eigenvalues of the operator dl 31 in the limit A, q — > and e n (k) are the corresponding 
eigenfunctions. The latter satisfy the orthogonality relation 



e n (k)D^ q (k)e m (k) = S nm (15) 

" A W k 



in the limit of small q and A. Using Eqs. (1 14-i and (1151 . Eq. (II Oi reduces to a simple geometric series and 
is easily summed, giving 

rf£ s (k,k') = i x X l PP( , e„(k)e„(k'). (16) 

n 1 - w 

Clearly, positive eigenvalues produce a singularity in the vertex when A is lowered to the energy scale A c 
given by A„B^ P (0) = 1, or equivalently 



2 See Section 5.7 of Nozieres' book 1301 for a comprehensive explanation. Rigorous mathematical statements about the moving 
Fermi surface have been presented in I31II32||3*3*1 . For recent investigations of interaction-induced Fermi surface deformations, 

To be more precise, the operator defined in Eq. 1131 is Hermitian with respect to the scalar product given in Eq. Q5). 
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The divergence of the vertex T BCS at A = A c is usually interpreted as signaling the onset of supercon- 
ductivity. In fact, the divergence is caused by low energy p-p pairs and requires some attractive channel 
in the interaction (corresponding to a positive eigenvalue of I BCS ). This suggests that the instability is 
associated with the formation of Cooper pairs. Moreover formula d!7l > for the critical energy scale corre- 
sponds exactly to the BCS formula for the critical temperature. To establish the connection to BCS-theory, 
consider the standard gap equation 




where g BCS (k, k') is the coupling function, which describes the scattering of Cooper pairs from mo- 
menta k and k towards k' and — k'. Close to the critical temperature, the mean-field dispersion E^ = 

A£ + £k is replaced by |£k|- The resulting equation coincides exactly with the eigenvalue equation of 
the operator {13) at q = and A = A c , if we identify -g BCS (k, k') with I BCS (k, k') and A k , up 
to a normalizing factor, with e„(k). We conclude that the eigenfunction e„(k) with the largest positive 
eigenvalue determines the momentum-dependence (s-wave, p-wave, etc.) of the superconducting gap. 

It is actually not necessary to calculate the complete function I BCS (k, k'), but only its value on the 
Fermi surface. In fact, the eigenvalue equation reduces at small energy scales to 

Xn e " (k) = W(o) £ S ^ lBCS ^ k ') e «( k ') + (19) 

Thus, the eigenvalues A„ and e n (k)\^ =Q are determined by the value of I BCS (k, k') on the Fermi surface 
only. In this sense, the dependence of I BCS on the energies £k and £k' can be classified as irrelevant. 
The same is true for the dependence on the frequencies fco and k' , which has been neglected. Only the 
dependence on the momenta moving along the Fermi surface is crucial. 



4 Superconductivity from repulsive interactions 

As an example, we will now calculate the superconducting instabilities for an extended 2D Hubbard model 
away from half-filling 4 . The electron hopping is restricted to nearest neighbors, i.e. £k = — 2t(cosk x + 
cos k y ) — /i. Fig.|4]shows the Fermi surface for a typical electron density. The variable 9 used to parame- 
terize the Fermi surface is defined as the radial angle. The interaction 

Hi = U rif r nj r + V n r n r i (20) 

r (r,r') 

contains a nearest neighbor term V in addition to the Hubbard U. This corresponds to a coupling function 
g(k 1 ,k 2 , k 3 ) = U + 2V[cos(k 3yX - fc 2 ,x) + cos(fc 3lJ/ - k 2 , y )]. 

We study the superconducting instabilities in two steps. First we calculate within perturbation theory 
the p-p irreducible vertex 1(9, 9 1 ) — / < ^5 s A=0 (k(6'), k(0')), where k(#) is the Fermi momentum corre- 
sponding to the angle 9. A discretization of the 9 variable into 32 patches is introduced to allow for a 
numerical treatment. The second step is the solution of the eigenvalue equation ( 1191 . This is a trivial 
numerical operation, once the (9-variable has been discretized. 

The eigenfunctions e n {9) can be characterized by their transformation properties with respect to the 
symmetry group of the square lattice. The D4 point group has five irreducible representations: A\, A2, B\, 
B2 and E. They are illustrated in Tabled 



Similar calculations were performed before by Zanchi and Schulz [ 38 | and more recently by Guinea et al. 1391 . Our results are 
in perfect agreement with theirs. 



Fig. 4 The Fermi surface of the nearest-neighbor tight-binding band for /i = —0.8 (i.e. n ft! 0.7). 



Irreducible representation 


Basis 


Simple cases 


Aj. 


cos Am9 


s-wave 


A 2 


sin Am9 




Si 


cos (Am + 2)9 


d*x 2 — y 2 


B 2 


sin (Am + 2)9 


dxy 


E 


( cos (2m + 1)9 
{ sin (2m + 1)9 


p-wave, /-wave 



Table 1 Irreducible representations of the D4 point group, a basis of functions, which span the corresponding sub- 
space (m is an arbitrary integer) and the term commonly assigned to the simplest order parameters in each case. The 
A2 subspace consists only of rapidly oscillating functions with at least eight nodes along the Fermi surface. 



v(0) X 

2.5 

2 

1.5 

1 ■ 
0.5 





Al 


* 


A2 


- ■ 


Bl 




B2 




E 



Fig. 5 The positive eigenvalues of 1(8, 6') multiplied by the density of states at the Fermi level as a function of the 
electron density for an attractive nearest-neighbor interaction V = — O.lt and U = t, as obtained by a first order 
calculation. The only positive eigenvalues are in the E and B\ symmetry block. The corresponding eigenfunctions 
are p-wave like (for E) and d a .2_j / 2-wave like (for B\). The quantity v(Q)\ determines the energy scale of the 
superconducting instability (see Eq. <17t ). 



All the calculations have been done for a repulsive on-site interaction U = t. In the case of an attractive 
V < 0, the appearance of superconductivity is not surprising. In fact, the first order vertex I = —g has 
positive eigenvalues in this case. The results can be seen in Fig. [5] They show triplet p-wave supercon- 
ductivity for electron densities n < 0.65 and singlet d^^-wave superconductivity for 0.65 < n (close 
to half filling, additional instabilities arise in the p-h channel). The second order contributions to 1(9, 9') 
give only minor changes to the present result, although they introduce small positive eigenvalues in the 
symmetry blocks B 2 ,Ai and A 2 . 



10 



B. Binz, D. Baeriswyl, and B. Doucot: Weakly interacting electrons and the renormalization group 



v(0 ^,v,0H 



0.06 



. 04 



. 02 






- Al 


- * - 


A2 


- ■ 


Bl 




- B2 




E 



0.2 



0.4 



. 6 



. 8 



n 



Fig. 6 Same as Fig.[5]from a second order calculation for the pure Hubbard model V = and U = 0.5t (Kohn- 
Luttinger superconductivity). The largest eigenvalue of each symmetry block is shown. The corresponding eigenfunc- 
tions are d^-wave like for B2 (n < 0.5) and d x 2_ y 2-wave like for Bi (n > 0.6), and approximately /-wave for E 
(which is leading only in a narrow range of dopings). The inset shows the global behavior with its drastic increase of 



the cL 



-wave eigenvalue close to half-filling. 



In contrast, in the pure Hubbard model (V = 0, U > 0), positive eigenvalues of 1(9, 9') appear only 
via second order terms. Since the sum of the three PH2 diagrams is zero in the Hubbard model, the only 
second order contribution to 1(9, 9') comes from the PHI diagram. The latter is computed numerically, 
using a mesh of 3000 x 3000 k-points in the Brillouin zone and a temperature T = O.OOli. Kohn and 
Luttinger have investigated the same diagram in their historical paper |40|. The difference is that we are 
dealing with an anisotropic 2D system as opposed to the isotropic 3D case considered in |40|. 

The results in Fig.|S]for the Hubbard model show d xy -wave symmetry for n < 0.5 and d 2 ,2„ J( 2-wave for 
n > 0.6. In between there is a narrow range of densities where triplet superconductivity (with an /-wave 
like order parameter) is favored. Note that the eigenvalues of Fig. [6] created by second order terms are 
in general smaller than those in Fig [5] with an attractive V, leading to superconductivity at lower energy 
scales. 

The appearance of superconductivity by second-order corrections to the vertex has a physical interpre- 
tation in terms of an effective attractive interaction, which is mediated by spin fluctuations. In fact, 1(9, 9') 
equals 

1(9,9') = -U-U 2 X o(k e +k e ,), (21) 

where xo(q) is me spin susceptibility of the non-interacting system. The attractive part of 1(9, 9') comes 
therefore from the spin susceptibility Xo(q)- It is thus fair to say that an effective attraction is created 
by the exchange of spin fluctuations. The situation is analogous to the conventional superconductivity 
of metals, where an effective attraction is created through the exchange of phonons, however with the 
important difference that the same electrons which feel the effective interaction are also responsible for the 
spin fluctuations. 

The idea of spin-fluctuation-induced pairing has led to semi-phenomenological theories of supercon- 
ducting materials with strong magnetic correlations such as the cuprates which are antiferromagnets in the 
absence of doping BT1l42ll43ll44ll45l. 
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These results strongly support that there is superconductivity in the repulsive Hubbard model at weak 
coupling. While for electron densities smaller than 0.6 we expect superconductivity only at exponentially 
small temperatures, the energy scale for d x 2_ y 2- pairing increases drastically as the system approaches 
half-filling. This effect is partly due to the enhanced density of states, but approximate nesting, which 
enhances spin fluctuations and thus the effective attraction among electrons, is likely to be more important. 
Very close to half filling, spin fluctuations become so strong that the ladder approximation is no longer 
valid, and a consistent treatment of both p-p and p-h diagrams is required. This is the subject of the 
remaining part of this paper. 

5 Nesting 

In Section [5] it was assumed that the divergences in the perturbation series of the vertex are exclusively 
generated by the p-p diagram. The divergence in the p-p diagram at k = depends on the parity relation 
£ p = £_ p , which is applicable quite generally. 

In contrast, divergences in the p-h diagrams arise only if the Fermi surface has the so-called nesting 
property. A Fermi surface FS is called p-h nested if the following conditions are satisfied. 

1 . There exists a vector Q such that FS has some finite overlap 5 with its own translation Q + FS. 

2. The overlap is such that the occupied part (£ p < 0) in the vicinity of FS is covering the empty region 
(£ p > 0) of Q + FS and vice versa. 

In this case the p-h diagram PHI shows a logarithmic divergence for qi = Q and PH2 is divergent for 
q2 = Q. For example every bipartite lattice (i.e. a lattice with electrons hopping only from one sub-lattice 
A to another sub-lattice B) leads to a p-h nested Fermi surface at half filling. 

Similarly one can define p-p nesting as follows. 

1 . There exists a vector K such that FS has some finite overlap with K — FS. 

2. The overlap is such that the occupied (empty) side of FS is covering the occupied (empty) side of 
K IS. 

In this case, the p-p diagram PP shows a logarithmic divergence for k = K. For example, every Fermi 
surface with spatial inversion symmetry is p-p nested for K = 0. 

These are necessary conditions for the occurrence of infrared divergences if the density of states at the 
Fermi level is finite, i.e. in the absence of Van Hove singularities. 

The situation is particularly rich if the Fermi surface has flat portions. In this case there is usually a 
continuum of nesting vectors satisfying the above conditions. RG or parquet calculations for this case are 
provided in Refs. 1461 1471 l48l . 

6 Bethe-Salpeter equations 

If the divergences arise in both p-p and p-h bubbles, the ladder approximation used in Section [5] is not 
sufficient. In fact, in such cases the leading term of the perturbative expansion (|8j is given by the so-called 
parquet diagrams. These diagrams are obtained by retaining the five one-loop (or second order) diagrams 
shown in Fig.H an d by replacing any bare vertex by one of the one-loop diagrams and continuing this 
process to any order. An example is given in AppendixlBl 



By "finite overlap", we mean an overlap with a finite (D — l)-dimensional volume. Thus a single overlapping point is sufficient 
for D = 1 (and ladder systems), for D = 2 it should be a curve of finite length, etc. 
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Fig. 7 The structure of diagrams, which are reducible in the p-h 1 channel (a) and the p-h 2 channel (b). The grey 
circles represent any sub-diagram. 

To keep track of all the logarithmic divergences, one defines three different "channels" of two-particle 
reducible diagrams. A diagram is called reducible in the p-p channel if it has the structure shown in Fig. [2] 
Similarly, the diagrams of the form shown in Fig. 0a and b are called reducible in the channels p-h 1 and 
p-h 2, respectively. 

Let us denote by 

R A p (ki, k 2 , ks) the set of reducible diagrams in the p-p channel, with an infrared cutoff A, 
R p A hl (hi , k 2 , k§) the set of reducible diagrams in the p-h 1 channel, 
R p ^ 2 [k\, k2,k^) the set of reducible diagrams in the p-h 2 channel, 
I\(k 1 ,k 2 ,k 3 ) the set of two-particle irreducible diagrams. 

It is rather simple to check that a given diagram is either two-particle irreducible or reducible in only 
one of the three possible channels p-p, p-h 1 and p-h 2. Hence 

T A = I A + Rf + Rf 1 + Rf 2 . (22) 

The set of irreducible graphs in each of the three channels o — pp, phi or ph2 is defined by 1^ = T A — R A - 
We already encountered I pp in Section|3| where it could be calculated within perturbation theory. This is 
no longer true in the present case, where infrared divergences occur in the p-h channels as well. 

Clearly, R A P is given by the ladder series Eq. dlOi . but without the first (zero bubble) term. The first- 
order term is missing since it corresponds to an irreducible contribution. We write this series formally as 



Rf = I PP D PP I PP + I PP D PP I PP D PP I PP + ... . (23) 
It satisfies the integral equation 

Rf = I PP D PP I PP + I PP D PP R PP (24) 
or, using the relation T = I pp + R pp , 

R pp = I pp D pp T K . (25) 
This is the Bethe-Salpeter equation 6 . With all the functional dependences included it reads 

R p A p (k 1 ,k 2 , k 3 ) = -^J2 J A P ( fc i. fc 2, k - p)D pp k {p) r A (p, k - P , fc 3 ), (26) 
p 

where k — ki + k 2 and Df k (p) = G A (p)GA(k - p). 




Fig. 9 The general structure of a reducible diagram in the p-h 2 channel. The wavy lines drawn horizontally stand for 
— 7^ h2 , but the wavy lines drawn vertically stand for — J? \ The electron lines stand for full propagators Ga. 



There are similar equations for the two other channels. First, R phl is given by a series over the p- 
h ladders shown in Fig. (HI where the wavy lines stand for —I A hl and the electron lines denote the full 
propagators Ga. In complete analogy with the p-p case one finds the Bethe-Salpeter equation for the p-h 1 
channel 

R A hl (k 1 ,k 2 ,k 3 ) = — ^I A hl (k 1 ,p + q 1 ,k 3 )D A h qi (p)T A (p,k 2 ,p + q 1 ), (27) 
p 



where q\ — k 3 — k\ and 

D p A h q (p)=G A (p)G A (p + q). (28) 

The p-h 2 channel is more involved. In fact, a general i? ph2 -diagram has the structure shown in Fig. |9j 
where the wavy lines drawn horizontally stand for —I A h2 , but the wavy lines drawn vertically correspond 
to —I^ 11 . The electron lines stand for full propagators Ga- Note that the diagram must have at least one 
horizontal wavy line (in the opposite case it is a p-h ladder and not reducible in the p-h 2 channel). 

The simplest examples of p-h 2-reducible diagrams are shown in Fig. [2 The first PH2 diagram of this 
figure is given by 

" 2 jv £ I ? 2 V ei >P>P + q ^°tm (P) + 92,k2,h), (29) 

v 



See for example I49II30I . The Bethe-Salpeter equation was originally introduced to calculate two-particle bound-states within 
quantum electrodynamics 1501 . 
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where q2 = — &2. We will write this specific convolution of the functions I p ^ l2 and as a formal 
product 

-2lf 2 Dflf 2 . (30) 

Note that this notation differs from the formal product introduced before in the p-p case. With this notation 
one can write 

RP h2 = _ 2 IP h2 D ph IP h2 + ( XI phi - )D ph I ph2 +I ph2 DP h( X phis j + " higher order terms", (31) 

where the overall A-index has been omitted for the simplicity of notation and XF(k\ 1 k2,k%) — F(k2, ^1,^3) 
for any function of three energy-momenta. The p-h ladder series can also be written using the same nota- 
tion: 

XR phl = (XI phl )D ph (XI phl ) + "higher order terms". (32) 
The following definitions turn out to be very useful. 

RC = 2R ph2 -XR ph \ 

I c = 2I ph2 -XI ph \ (33) 
T c = I c + R c = (2- X)T, 

and 

R s = -XR ph \ 

I s = -XI ph \ (34) 
V s = I s + R s = -XT. 

The superscripts c and s refer to charge and spin, respectively. Note that T c / S = ± T^, so these vertex 
functions enter naturally in the calculation of charge- and spin- response functions. 
Equations ( 13 11 and d32l > can now be written in the simple form 

RS ,c = _j8, CD phjs,c + "higher order terms". (35) 
The exact expressions for R s,c are given by the infinite series 

R s,c _ _js,cjjphjs,c ^ js,c jyph js,c jyph js,c 

In order to prove Eq. d36i . we first write R ph2 as an infinite series R ph2 = 52?Li R p I l2 , where R P t h2 is the 
set of diagrams, which are exactly ?i-times reducible in the p-h 2 channel. A similar decomposition is done 
for R phl . The recursion relations between , R P ^ 2 X and R phl , R ph2 are shown graphically in Figs. 1101 

and^2 The relation for RfJ 12 CFig. 1101 can be understood as follows. A general .^-diagram as shown 
in Fig.|9lcan either start on the left with a vertical or a horizontal wavy line. If it starts with a vertical line, 
then it is given by the second diagram of Fig. If it starts with a horizontal line, there are two cases. 
If there are no other horizontal lines except for the starting one, the diagram is given by the third term of 
Fig- El If there is more than one horizontal line in the diagram, it is given by the first diagram in Fig. 1101 
Analytically the equations depicted in Figs. 1101 andll llread 

flpM _ _ 2I ph2 D ph R ph2 +XI phl D ph R ph2 +I ph2 D ph XR phl (3?) 

XR P ^ = XI phl D ph R phl , (38) 
or, with the definitions J33I and J34> . 

R^ = -I sc D ph R s n c . (39) 

From this, it is easy to deduce Eq. i36\ . In analogy to the p-p case, we find now the Bethe-Salpeter 
equations for V s and T c 

R s,c = _ IS ,c D ph T s,c_ (40) 
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Fig. 10 The recursion relation for i?^, 2 j . Wavy lines drawn horizontally stand for I ph2 and wavy lines drawn verti- 
cally stand for I p . The grey rectangles with index n stand for R„ h2 (R^) if they are drawn horizontally (vertically). 



n+l 



k, 



Fig. 11 The recursion relation for . The conventions are as in Fig. llQI 



7 RG flow of the vertex 



The Bethe-Salpeter equations J26b and J40i . three integral equations for the three unknown functions R pp , 
R° and R s , are in general difficult to solve. We will address the more modest task of calculating the 
leading term in the perturbative expansion l|8). This will be accomplished via the flow equation, i.e. a 
differential equation for 9aT where we keep only the leading terms in A — > 0. Within this approximation, 
the two-particle irreducible vertex / is given by the bare interaction — g. Thus by Eq. J22L 



r (h , k 2 , k 3 ) = Rpp (h , k 2 , fc 3 ) + R* hl (h , k 2 , fc 3 ) + R ph2 (fci , fe, fc 3 ) , 



(41) 



where the dot means a partial derivative with respect to A. 

Consider the first term R pp . The derivative of the Bethe-Salpeter equation d25l > has three contributions 



R pp — jppjjppy -f pP jjppy + JPpjjppy 



(42) 



Although I pp is by no means negligible by itself (since it contains terms which are reducible in the p-h 
channels), the first contribution, IPPD PP T, can be shown to be of subleading order in A. We will not prove 
this here, but an example is discussed in detail in Appendix 151 The last term is written as Ippj^ppjpp + 
Ipp DPP rpp an d iPPD pp iP p is neglected for the same reason as i pp D pp T. Therefore 



RPP = jpp jjppy 4 jpp]jppjipp_ 
This equation can be iterated to give 

rpp = jppjjppy 4 i p p jypp i pp d pp t 4 pp jypp pp jypp pp jyp p Y - 

= JPpjjppy 4 Rppjjppy 
= TDP p T, 

or written out with the full functional dependencies 

R p 1 P(k 1 ,k 2: h) ^ — ^T A (k u k 2: k- p)D^ k (p) T A (p, k-p,k 3 



(43) 



(44) 



(45) 



where k — k\ + k 2 is the total frequency-momentum. 
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c* 



it 2 

jphl 




Fig. 12 Negligible terms in the flow equation for <9aF. The grey circles stand for any sub-diagram. 

The same kind of differential equations are obtained in an analogous way for the p-h channels 7 . Using 
Eqs. and 

we obtain 

RS,c = _ I s,cf ) ph T s,c _ p.cjjphfts^ (46) 

js ,c jjph-^s ,c _|_ js ,c jjph js ,c j^ph-ps ,c (47) 

= _J8,CQphY8,c _ j£S,C£)ph-ps,c ^4g-j 

_ _Y s ' c £)P fl Y SC (49) 

The only point to be verified is that the terms j s > c jjP h Y s > c and I s ' c D ph I s ' c , which have been left out in 
Eq. j46L are in fact of subleading order in A. Writing out I s c in terms of I phl and I'P h2 , one can verify that 
all the neglected terms are of the form shown graphically in Fig.^1 They can be shown to be negligible 
in analogy with the example of AppendixIBl Equation i49i , expressed in terms of R phl and R ph2 leads to 

R phl = X ((XT)D ph (XT)^J (50) 
R ph2 = ^TD^T+iX^D^T + TD^iXT), (51) 
or written out: 

R P ! ll (k 1 ,k 2 ,k 3 ) = -^^2D P ^ qi (p)r A (p,k 2 ,p + qi)T A (k 1 ,p + q l7 k 3 ) (52) 

v 

Rf 2 (k u k 2 ,k 3 ) = -^^D p A h ; q2 ( P )[-2r A (k 1 , P , P + q 2 )r A {p + q 2 ,k 2> k 3 ) (53) 

p 

+T A (p, ki,p + q 2 )T A (p + q 2l k 2 , k 3 ) + T A (ki 1 p,p + q 2 )T A (k 2 ,p + q 2l k 3 )} , 

where qi = k 3 — k x , q 2 = k 3 — k 2 are the direct and exchanged transfered momenta and D A ' l (p, q) = 
G A {p)G A {p + q). 

Equations J41i . ( 145 > . \52\ and (1531 are the one-loop RG equations we were looking for. They describe 
the behavior of the vertex under a differential change of the energy scale A. For the sake of clarity, these 
equations are collected as follows: 

f (kx, k 2 ,k 3 ) = R pp (k u k 2 , k 3 ) + R phl {k 1: k 2 ,k 3 ) + R ph2 (k 1 ,k 2 , ha) (54) 

R p A {ki,k 2l k 3 ) = — ^2 D p A k{p) r A (fci, k 2 ,k- p)r A (p, k-p, k 3 ) 
p 

R P A hl (ki, k 2 ,k 3 ) = — ^ D p A h qi (p) T A (p, k 2 ,p + q 1 )T A (k 1 ,p + q x , k 3 ) 



1 Remember that the formal products of functions mean something different in Eq. 1491 than in Eq. 1441 . The formal product in 
the p-p channel is defined by Eqs. 1251 and 1261 . whereas the formal product in the p-h channels is defined by Eqs. 1291 and 

OoJ. 
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R P A h2 (k u k 2 ,k 3 ) = —Y,D p A h q2 (p)[-2r A (k 1 ,p,p + q 2 )r A (p + q 2 ,k 2 ,k 3 ) 

p 

+T A (p, ki,p + q 2 )T A (p + q 2 ,k 2 , k 3 ) + T A (ki,p,p + q 2 )T A (k 2 ,p + q 2l k 3 )} , 

where k = k\ + k 2 is the total frequency-momentum and q\ = k 3 — k\, q 2 — k 3 — k 2 are the direct 
and exchanged transfered frequency-momenta, respectively. The p-p and p-h propagators are defined as 
D P A P k (p) = G A ( P )G A (k - p) and Df q {p) = G A (p)G A (p + q). 

These equations have a very simple diagrammatic interpretation. The three terms R pp , R phl and R ph2 
correspond to the one-loop diagrams shown in Fig.^ where the wavy lines now represent full vertices T A 
and the product of two single particle propagators has to be derived with respect to A. Since D PP ^ ph (p) = 
G A (p)G A (q ± p) + G A (p)G A (q ± p), each diagram can be viewed as the sum of two terms, where one 
of the two lines represents the propagator G A (p) and the other is G A (p), the propagator restricted to the 
energy scale A. 

Although we have defined A as a sharp infrared cutoff in the band energy, the RG equations are also 
correct for a cutoff in the frequency or for smooth cutoffs. It is straightforward to use the finite temperature 
instead of a cutoff to regularize the theory. The RG equations for Ot^t are then obtained by replacing in 
Eq. (HI 

"4 £ d T (±D«(p,q)\ po=2 „ n2n+1) )..., (55) 
^ p neZ, p ^ ' 

where o = pp, ph. 

Different regularizations of the theory give the same result to leading logarithmic order. For instance, 
the vertex at zero temperature and finite infrared cutoff A is, within logarithmic precision, equal to the 
vertex at temperature T = A without cutoff. 

Let us close this section with the remark that the present RG formalism allows the calculation, to leading 
order in the infrared cutoff, of various generalized susceptibilities, i.e. the linear response of the system to 
weak external perturbations. The required relations are derived in Appendixlcl 



8 Discussion of the low-energy behavior 

The leading terms for A — * of the right-hand side of the one-loop RG equations behave like ~ r A A _1 
(note that A -1 is the derivative of the logarithm). However, the one-loop RG equations contain in general 
many non-leading terms, like R A P if k is not a p-p nesting vector and R ? 1 if q! is not a p-h nesting vector, 
etc. These terms are diverging at most as ~ T\ log A. 

Since the one-loop RG equations are only correct to leading order in A, only the leading terms are to 
be taken seriously. Taking into account subleading terms in the one-loop RG equations is inconsistent and 
arbitrary, since these subleading terms are comparable to terms that have been neglected. The contributions 
to the exact RG flow can be classified into three classes: 1. leading terms of the one-loop equation, 2. 
subleading terms of the one-loop equation and 3. terms which are neglected in the one-loop equation. One 
has to compare the subleading terms ~ T\ logA with the neglected terms. According to a careful analysis 
by Salmhofer and Honerkamp |51 1, the first neglected terms are or the order ~ T\ log 2 A (see Eq. 108 of 
ED)- Schematically, 

Fa = "leading terms" + "subleading terms" + "neglected terms" . (56) 

V v ' V v ' S v ' 

~r|A-! ~r2i og A ~r3io g 2 A 

The RG flow is split into three different energy regimes. 
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1. High energies, in which A is not small. In this regime, the subleading terms of the one-loop equation 
dominate the neglected two-loop terms, provided Ta logA is small. Thus the one-loop RG equations 
can in principle be used to calculate T\, starting with a cutoff equal to the bandwidth, where the 
vertex is given by the bare coupling. In this regime, there is no small parameter to further simplify the 
functional RG equations, so using them is technically very difficult. But in the high energy regime 
(g log A <C 1) naive perturbation theory provides a controlled and much more feasible method of 
calculation than the RG. 

2. Small energies, where A is small such that (A log A) -1 3> Ta logA ~ 1. Naive perturbation theory 
breaks down in this regime, but the leading terms of the one-loop RG equations are superior to the 
neglected terms even if T is not necessarily small. This is where the one-loop RG is most useful. 
Note, that the subleading terms are no longer superior to the neglected terms in this regime. Thus, 
they have also to be neglected to be consistent. The small parameter A justifies to replace every vertex 
T(ki, k2, fe) in the one-loop RG equations by its value on the Fermi surface (fe u ,i = &o,2 = fco,3 = 
£ki = £k 2 = £k 3 = 0), for reasons given in the end of Section[3](for the case where the instabilities 
only occur in the p-p channel) and in AppendixlBlffor the nested case). 

3. The critical regime, where A is close to the critical energy scale A c at which T\ is diverging. If A 
is too close to A c , the neglected terms are no longer negligible and the one-loop approximation is no 
longer accurate. Note however that the weaker the initial interaction is, the more A can approach A c 
before the one-loop approximation breaks down. 

This remark about subleading terms concerns in particular the self-energy corrections to the propagator 
Ga- As already pointed out after Eq. dlOi . the self-energy changes the shape and location of the Fermi 
surface. This effect is by no means negligible in perturbation theory. Other effects, such as the renormal- 
ization of the Fermi velocity or the reduction of the quasi-particle weight are of subleading order. So these 
effects cannot be taken into account consistently within the one-loop approximation. For this reason, it is 
a consistent approximation to replace the dressed electron propagators Ga in the one-loop RG equations 
by the bare ones C\. 

9 Relation to the Wilsonian approach 

A key ingredient to the Wilsonian RG is the idea to replace the given problem by a different one with less 
degrees of freedom but with the same low energy behavior. To achieve this, the effect of the eliminated 
degrees of freedom is incorporated in a renormalization of the parameters of the effective low-energy 
theory. This strategy has been successfully followed using Brillouin-Wigner perturbation theory in the 
strong coupling limit of various many -body problems (see for example 1521 1531 l54l ~). The best known 
example is the Hubbard model at half-filling which is represented, in the limit of strong coupling, by a 
Heisenberg spin Hamiltonian in the limit of strong coupling. In the weak coupling limit, the application 
of the Brillouin-Wigner formalism to compute an effective Hamiltonian is less evident and has not been 
followed to our knowledge. A tractable implementation of the same idea (i.e. elimination of high energy 
degrees of freedom and renormalization of parameters of the effective theory) is given more easily in the 
functional integral representation 8 . Some basic definitions and results of the functional integral formalism 
are given in ApDendixIDl for a detailed presentation, see |60|. 



It should nevertheless be mentioned that the idea of effective Hamiltonians on a reduced Hilbert space led to a most powerful 
numerical tool for one-dimensional systems: the density matrix renormalization group (see [55] for an introduction). An 
alternative route to effective Hamiltonians was presented recently by Wegner 1 56 ] and applied to the 2D Hubbard model 1571 
and other correlated systems 1 581 1591 . 
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The bare propagator is endowed with an infrared cutoff A on the band energy C\(k) = 0(|£k| — 
A)C(fc). The effective interaction depending on the energy scale A is defined as 

W A lx\ = - log / d^c A bP} e- w ^ , (57) 



where d^c A M * s me normalized Gaussian measure according to the propagator C\ (Eq. dl30» . and W 
is the model interaction in terms of Grassmann variables Eq. J124i . Note that the integration with respect 
to dfic A bP] i s perfectly defined, although C^ 1 is not. This can be seen most easily in the expansion of 
Wa[x] in terms of Feynman diagrams. The evaluation of these diagrams involves only C\ and never C A . 
Whenever appears in an intermediate step of a calculation, it may be regularized by replacing the zero 
in the Heavyside function by an infinitesimal number. 

The effective interaction Wa [x]> which depends on the Grassmann field x, has a twofold interpretation. 
On the one hand, we can restrict the field \ to the low energy degrees of freedom ip^ = 0(A — l^kDV'fco-- 
The object 

sfbP<] = crV) - w A [^<] (58) 

corresponds then to Wilson's effective action, which describes the system in terms of i/j < only. In fact, for 
observables (or Green's functions) which depend only on the low energy fields, one can show 

(0[^<}) = \\ ViP< j Vip> e s[ ^ < ^ >] 0[ip<] (59) 

1 

Z 

and 

Z= [v^< [ V^> e 3 ^^ = I V^< e 3 "^. (61 ) 



\ I V^< e 3 °"^ 0[^<] (60) 



On the other hand, Wa is the generating functional of amputated connected correlation functions with 
infrared cutoff A because of the identity 

log Z K [ v ] = -(fj,C A r])-W A [C A rj\ , (62) 

where Z A is given by Eq. dl29> . with C replaced by C A . To prove Eq. J62l >. we write out the Gaussian 
measure in (Eq. dl30» explicitly. Then, only in the numerator, we perform a shift of the integration 
variables ip — ► i/j — \- Equation J62l > follows finally by putting \ = G A r\ and taking the logarithm. 
As a consequence of Eq. i62\ . the quadratic part of Wa is related to the self-energy Ea by 

, S ,. WaM = -C^(k)-(yj ak ^ k ) A C A 2 (k) (63) 

i-c A (fc)s A (fc)' K } 

where we have used the following identity for the full electron propagator G A (k) = — (tpak^ak) A = 
C A (k)(l — C A (k)'S A (k))~ 1 . Therefore in the case |£k| < A the right-hand side of Eq. i64i simply 
becomes T, A (k). 

Similarly the quartic part of Wa is related to the one particle irreducible vertex Ta- In fact, differenti- 
ating Eq. (I62> we find 

4 

-(^ak 1 4'a'k 2 'lpa'k 3 , ipcrk i )c,AY]_ C A 1 ( k i) ( 65 ) 



5 4 W A [x] 



SX<jk i SXv'k 3 SX>j'k 2 SXcrk 1 



x=o i=l 

rr'(fci,...,fc 4 ) 

^ni=i[ i -^A(Ai)SA(*i)]' 
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In the last line we have used Eq. dl27> and the definition of the vertex (Eq. ( I78t ). 

The quartic part of Wa is of the same form as Eq. (I124> with an effective coupling function g&(ki, k2,k 3 ) 
that now depends on the frequencies as well as the momenta. Taking functional derivatives of Eq. J124i . 
we find for |£ kt | < A, 

(1 - 5 aa ,X)g K {kxM, h) = -IT'(fci, k 2l fes) (67) 

and thus 

g\(ki 1 k 2 ,k 3 ) = -T A {ki,k 2 ,k 3 ). (68) 

We summarize that g\ is equal, up to the sign, to a connected amputated correlation function if all |£k, | > A 
and to the one-particle irreducible vertex in the opposite case. g\ is therefore not continuous at |£k ; | = A. 
A formal and non-perturbative proof of these relations was given by Morris 1 6 1 1 for a bosonic field theory. 
The derivation given above is perturbative, but a generalization of the non-perturbative proof of Morris to 
fermions appears to be straightforward. Morris has also shown that Ea and Ta are continuous at |£k 4 1 = A, 
in contrast to g\. 



10 Other one-loop RG equations 

The effective interaction satisfies the following exact RG equation for d\ Wa = Wa: 

*-t OXakOXak 'H OXak kick 

a,k (T,k 

This equation was first derived by Polchinski in the context of a scalar field theory 1 62 1 . Zanchi and Schulz 
1241 125 1 proposed to develop Wa up to order six in the fermionic variables and to neglect terms of higher 
order. Terms of order six are not present in the original interaction but they are produced by the RG 
procedure. Their effect is then to renormalize the effective coupling function g\. The result is a closed 
one-loop equation for the coupling function g\(ki, . . . , k^) where |£k< | < A. It is identical, within the 
correspondence g\ = — T\, to our Eq. i54i with one difference. In the RG equation of Zanchi and Schulz, 
the vertices on the right hand-side are not evaluated at the scale A, but at a higher scale A, which is given by 
the band energy of the single-particle propagators, i.e. A = Max{|£ p |, |£k- P |} in the p-p term of Eq. J54> . 
A = Max{ |£ p | , |£ p + qi | } in the p-h 1 term and A = Max{ |£ p | , |£ p + q2 1 } in the p-h 2 term. Their equation 
is thus non-local in A. It is a flow equation with memory, i.e. the flow at scale A doesn't only depend on 
the vertex function Ta but also on the history of the flow. 

Since this is not very convenient, it was proposed 1 63 1 to develop Eq. j69t into Wick ordered polyno- 
mials of the fermionic variables instead of monomials as it was done above. Wick ordering with respect 
to the low energy propagator D\ = C — C\ results in the same one-loop equation as above but now all 
the couplings are evaluated at the actual RG variable A and — d [Cj^(p)C\(q)] /dA has to be replaced by 
d [Da(p)-Da(<3')] /dA, i.e., one propagator is at the energy A and the energy of the second propagator is 
now restricted to be smaller than A. This is different from our one-loop RG equation, where the second 
propagator is restricted to higher energies. 

We have seen above that the coupling function of the low-energy effective theory equals, up to a sign, the 
one-particle irreducible vertex T\. This correspondence can not be generalized to higher order vertices, i.e. 
it would be completely wrong to say that the sixth order term of Wa is related to the one-particle irreducible 
part of the three particle Green's function and so forth. In addition, the correspondence Eq. d68l relies on 
the choice of the cutoff which ensures that C\(p) = 0, for |£ p | < A. The correspondence no longer holds 
for an alternative scheme, where for example the finite temperature is used to regularize the theory instead 
of the infrared cutoff. 
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The general one-particle irreducible vertices are obtained by performing a Legendre transformation on 
the functional Wa- Wetterich has presented a renormalization group scheme for bosonic field theories, 
working with this Legendre transformed quantity rather that with the effective action defined in Eq. J57> 
1641 1651 . The idea was implemented recently for the many-fermion problem in 1271 1511 . The resulting 
RG equation is identical to Eq. (I54> apart from self-energy terms which should be, as we argued above, 
neglected for consistency reasons. In fact, if the cutoff is introduced by some function multiplying the bare 
propagator 

Ca(p) = t^-, (70) 

then the full propagator is given by 

G A (P) = : r CA ( iP l p ; (71) 
Oa(p) 

iPo - - ^a(p)Sa(p)' 
The derivative with respect to A gives 



(72) 



G A ( P )= "JM^M ■ °)^— 2 . (73) 

The first term of the right hand-side is the single-scale propagator of Ref . 1 27 1 . The one-loop RG equations 
of Salmhofer and Honerkamp are obtained from our Eq. J54l > by omitting the second term of Eq. (1731 . 

The RG equations of the three groups 1251 1631 127 1 differ in the treatment of the diagrams which give 
not a leading order contribution in the limit A — > 0. For example, the (leading order) p-p diagram at total 
momentum k = features two internal propagators with exactly opposite momenta. They have exactly the 
same band energy. So the non-locality of the Zanchi-Schulz equation is not present in this diagram. For 
the same reason, the energy-constraint of the Wick-ordered scheme can not have any effect in this case. 
So the three RG schemes treat such diagrams identically. The same is true for the p-h diagrams in the case 
of perfect nesting. However, as we have argued in Section only the leading order terms of these RG 
equations make sense. The subleading terms are of the same order as others that have been neglected. In 
fact, a consistent treatment of subleading terms requires going beyond the one-loop approximation. 

The equivalence of the different RG equations to leading order can also be understood in the following 
way. By successively integrating the one-loop RG equations (Eq. i54l ) and expressing the result in terms 
of Ta ~ —g, one obtains the full series of parquet diagrams. However the structure of these RG equations 
introduces a constraint on the energies of internal lines. For example the parquet diagram of Fig. 1131 is 
generated by the one-loop RG with the following constraint on the propagators 1,2,3 and 4 

Mm{|^|, |e 2 |} < Min{|e 3 |, l^l}- (74) 

Higher order parquet diagrams are generated with similar constraints, i.e. with a certain "ordering" of the 
band energies, when one goes from the inner loops to the exterior loops of a given parquet diagram. It can 
be checked by introducing this constraint into Eq. d86l > of AnnendixlBl that the constraint does not change 
the value of the diagram to leading logarithmic order in A. 

The different RG equations 1251 1531 1271 all generate the whole series of parquet diagrams, but the 
constraints are different. In the Wick ordered scheme, the constraint 1741 is changed into 

Max{|a|,|6|}<Max{|£3|,|£4|}. 

The RG equation of Zanchi and Schulz introduces the most restrictive constraint 

Max{|a|,|6|}<Min{|£ 3 |,N}, 
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Fig. 13 Fig. a) shows an example of a two loop parquet diagram, where the internal electron propagators are numbered 
from 1 to 4. Fig. b) shows the same diagram, but the wavy interaction lines are contracted to points in order to make 
the loop structure more transparent. 

i.e. both propagators of the inner loop are higher in energy than both propagators of the exterior loop. All 
these constraints are irrelevant for the leading logarithmic order in A. We conclude that our one-loop RG 
equations and those of Refs. 1 25 63 27 1 are all equivalent to leading order in A. 

11 Conclusion 

In conclusion, we have presented the RG method for weakly interacting electrons to leading (one-loop) or- 
der using a simple framework. We have chosen the point of view of partial summations of the perturbation 
series, taking into account the leading infrared-diverging terms and neglecting subleading contributions in 
an unbiased and consistent way. The procedure is controlled by two small parameters, the bare interaction 
strength and the energy scale of interest. 

In the case of a generic Fermi surface in more than one dimension, this procedure leads to the ladder 
approximation. Our calculations for the two-dimensional Hubbard model away from half filling show a 
pole in the scattering vertex T(pi,p2,P3) at pi + p2 = 0, where pi.2 are the frequency-momenta of the 
incoming particles. The standard interpretation of this pole is the development of a two-particle bound 
state, i.e. of a pairing instability. The symmetry of the pairing state is e?-wave. The pairing mechanism can 
be interpreted in terms of an effective attraction between electrons which is created by spin fluctuations. 
The critical energy scale at which the pole is created increases drastically as half-filling is approached. The 
reasons are twofold. First, antiferromagnetic spin fluctuations are enhanced due to approximate nesting 
and second, the density of states increases due to the proximity of a Van Hove singularity. 

However, if the Fermi surface is nested, the summation of p-p ladder diagrams is no longer sufficient, but 
the more general RG formalism involving p-h as well as p-p contributions is required. As opposed to some 
more modern formulations I TT61 1251 l63l 1271 . we have derived the one-loop RG equations by identifying 
the leading logarithmic terms. We feel that the less modern route chosen here makes the meaning of the 
one-loop approximation more transparent. In particular, it gives a natural guideline for deciding which 
terms are to be included in the calculation and which ones are to be neglected for reasons of consistency. 
For example, we argue in Section[8]that if the Fermi surface is not nested and the Fermi level is not close 
to a Van Hove singularity, the one-loop RG is not better than the ladder approximation. This is also true in 
the presence of Umklapp scattering, in contrast to the assertions of Ref . [ 27 1 . Furthermore, in Section fTol 
we compare three different versions of RG equations, which were discussed in the literature. We find that 
all of these approaches are equivalent within the regime of applicability of the one-loop approximation. 

The RG approach is most useful in one dimension, or in the case of a flat Fermi surface with many 
nesting vectors (see Section [5] for a precise definition). If there is only one single nesting vector, the 
situation is simpler. Consider for example an anisotropic square lattice at half filling. The dispersion is 
given by £k = — 2(t x cos k x + t y cos k y ), where < \t y /t x \ < 1 (not too close to neither or 1). Because 
the Fermi surface is perfectly nested by (it, 7t), there are infrared singularities in the p-h as well as in the 



23 



p-p channel. However, the non-ladder parquet diagrams appear not to be of the leading order in this case. 
As a consequence one can just sum the p-h 1, p-h 2 and p-p ladders separately. In other words, the one- 
loop RG equations decouple (to leading order) into three different sets, related to pairing-, spin- and charge 
instabilities, respectively |66 1. The study of the p-h instabilities is then in perfect analogy to our treatment 
of superconductivity in Section [5] and |4] A repulsive Hubbard-like interaction will lead to a spin-density 
wave instability. 

More difficult and also more interesting is the case, where the Fermi surface passes through saddle 
points of the dispersion relation. Because the dispersion is quadratic rather than linear at the saddle points, 
they render infrared divergences more dangerous, particularly in two dimensions, where saddle points lead 
to a logarithmic Van Hove singularity in the density of states. 

We have recently studied the half-filled nearest-neighbor tight-binding band on a square lattice (i.e. 
the case mentioned above, but with t x = t y ), where the Fermi surface is a perfect square and contains 
two saddle points. Our investigation yields a rich phase diagram as a function of the model interaction 
|67|. Apart from s- and d x 2_j,2-wave superconductivity, spin- and charge-density waves as well as two 
phases with circulating charge or spin currents (d-density waves) appear as possible instabilities. Unlike 
numerical studies of the functional RG equations 1251 163 , 68 1 which rely on a discretization, we do not 
obtain a simultaneous divergence of different, say superconducting and spin-density wave susceptibilities 
as the energy scale is lowered, but only of the leading susceptibility. Although there is no decoupling of 
the RG equations into independent sets, it turns out that the leading channel asymptotically decouples from 
the others as an instability is approached |67). 

The case of the half-filled nearest-neighbor tight-binding band is special because it features perfect 
nesting and saddle points at the same time. Considerable effort has been spent on the more general case of 
a Fermi surface which is at Van Hove filling but not nested. This case is particularly interesting because 
it potentially contains new physics such as non Fermi liquid behavior |69 70 1, itinerant ferromagnetism 
1711 1721 or Pomeranchuk instabilities |26|. Unfortunately, the one-loop RG equations are likely to be 
unreliable in this case. The reason is that the one-loop equation generates a correct summation merely 
of the leading divergent terms in the perturbation series, but the above-mentioned phenomena depend on 
subleading logarithmic divergences. This can be seen by considering the following example. To study the 
ferromagnetic (Stoner) instability, we have to calculate the scattering vertex T(ki, ka, k^) for an infinitely 
small momentum transfer, k/3 — ki — > 0. The instability is usually created by p-h ladder diagrams (Fig. [8}. 
In the case of Van Hove filling, the zero-momentum limit of the p-h bubble 9 diverges as ~ log T for low 
temperatures 10 as a consequence of the infinite density of states. This is a subleading divergence compared 
to the behavior of the zero-momentum p-p bubble, which is ~ log 2 T in the case of Van Hove filling. The 
third-order diagram of the p-h ladder series behaves as U 3 log 2 T since it contains two low-momentum 
p-h bubbles, each contributing with one logarithm. The one-loop RG theory in addition generates parquet 
diagrams, where p-p and p-h bubbles are mixed. The diagram shown in Fig. ^] contains one p-h bubble 
and one (internal) p-p bubble. If the external legs are chosen close to the saddle points, its leading behavior 
is U 3 log 3 T. The error from the one-loop approximation is of subleading logarithmic order, i.e. in our case 
of the order U 3 log 2 T. The problem is obviously that the error produced by the one-loop approximation 
is of the same order as the term which produces the ferromagnetic instability. The difficulty is not only 
encountered for third-order diagrams but proliferates to every order in perturbation theory. We conclude 
that the one-loop RG is an uncontrolled approximation for discussing ferromagnetism. 

A consistent weak-coupling treatment should at least include a self-consistent determination of the 
Fermi surface, self-energy corrections to both the quasi-particle weight and the Fermi velocity as well as 
two-loop corrections to the RG equations of the vertex. Although some progress has already been made on 



The p-h bubble is defined by B ph (q) = l/(/3V) ~^2 p C(p)C(p + q) and is equal up to the sign to the non-interacting spin 
susceptibility. 

' It is not possible to study the long-wavelength p-h response in the presence of an infrared momentum cutoff. The regularization 
has therefore to be done in another way, for example by a finite temperature. 
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certain of these aspects 1731 1741 l36l . a satisfactory control of the Van Hove case has not yet been obtained. 
New techniques will have to be developed to settle these issues. 
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A Definition and properties of the vertex function 

In this appendix, we give the definition of the imaginary-time two-particle Green's function and the vertex 
function. These quantities are not directly observable but they are analytically related to various physical 
response functions B49ll75ll60l . 

The two-particle Green's function is defined as 

G™' (n , r 2 , t 3 , ka , k 2 , k 3 ) = V < Tc CTkl ( n ) c CT , k2 (r 2 ) 4, k3 (r 3 ) 4 ki+ka _ k3 (0) ) (75) 

where V is the system volume (i.e. the number of lattice sites) and T is the time ordering operator. The 
imaginary time dependence of any operator is given by 

0(t) = e T{H -^ O e- T ^ H -^ , (76) 

where /i is the chemical potential and N — J2ka c k<r c ker tne particle number operator. The Fourier trans- 
form with respect to the r variables, 

rP 

G°° (h,k 2 ,k 3 ) = / dT 1 d75d7ae ,(fcoiri+fcoa7a - fcM ^G^ r (ri,75,7a,ki,k 2 ,k3), (77) 
Jo 

is a function of the D + 1-dimensional frequency-momentum vectors fe$ = (&rjj,kj). The Matsubara 
frequencies k^ run over the odd multiples ofn/P, because G aa is anti-periodic with period j3 in each of 
the r variables. 

The vertex function T aa (fci, k 2 , k 3 ), defined as the one-particle irreducible part of G aa (fci, k 2 , k$), 
satisfies the relation 

4 

G™'(k u k2,k 3 ) = f3V(6 k2M - d^5 kltk3 )G(k 1 )G(k 2 ) +r™'(k 1 ,k 2 ,k 3 )l[G{k j ), (78) 

i=i 

where G(k) is the exact one-particle Green's function and k^ = k\ + k 2 — k 3 . 

At this point it is worthwhile to investigate the symmetries of the vertex function. A simple SU (2) 
transformation shows that for a spin-rotation invariant system, 

T TT = (l-X)T n , (79) 

where XT(ki, k 2 , k 3 ) — T(k 2 , ki, k 3 ). Because the whole information is contained in the function T^, 
we will only consider the vertex of anti-parallel spins from now on and omit the spin indices (i.e. T = F**). 
The permutation symmetry gives 



F(kx, k 2 ,k 3 ) = F(k 2 , fa,ki + k 2 - k 3 ). 



(80) 
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q+Q 




k'-Q q 
Fig. 14 Example of a two loop parquet diagram. 



q-p+k 



k-Q 



Time reversal invariance, spin rotation invariance and parity imply 

T(ki, k 2l k 3 ) = T(fci + k 2 - k 3 , k 3 , k 2 ). 
From the behavior under complex conjugation one obtains 



T(ki,k 2 ,k 3 ) = T(fci + k 2 - k 3 , k 3 ,k 2 ), 
where k = (— ho, k). It follows that V is real if the frequencies are put to zero. 



(81) 



(82) 



B Parquet diagrams: An example 

A typical parquet diagram, is shown in Fig.^J It is a p-p diagram embedded into a p-h 2. For simplicity, we 
choose the total momentum and the external frequencies ko, k' and Qq to be zero and we assume a constant 
coupling function g. In addition, we assume that Q is a perfect nesting vector, such that £q+Q = — £ q . The 
value of the diagram (at T — > 0) is then 



2g 3 



(pvy 



J2 CK{q)C A {q + Q)C A (p)C A (q -p+k). 



After integration with respect to po and qo, we obtain 



©(fp£q-p+k) 



|£p| + |£q-p+k| 



w 



d£ 



2A 



KKKI + fc'l)' 



(83) 



(84) 



where the sum is over all momenta p and q allowed by the infrared cutoff and 



v 2 ^ 

pq 



s(t - - e P - eq-p + k)e(c P eq-p + k). 



(85) 



If the Fermi level is not at a Van Hove singularity, is usually bounded. Depending on the 

situation, v could even be suppressed to zero for £, £' — > 0. This would reduce the singularity of the 
diagram in such a way that it can be neglected within the leading order approximation. The cases where 
parquet diagrams are most singular are those where z>(£, £') is positive for £, £' — > (and can thus be 
approximated by a constant). It means that to leading order, our diagram is proportional to the following 
integral: 



P(A) 



w 



df 1 , 2 W 



(86) 
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For further convenience we have introduced the symbol Fa(£) for the value of the internal p-p bubble. The 
result shows that indeed our diagram is of the leading order g 3 log 2 A. 

To obtain its contribution to the RG equation, we take the derivative with respect to A 

-^(A) = ^a(A) + / a "^ (87) 

The first term is obtained by putting the large loop variable £ to the energy scale A. The second term comes 
from the derivative of the internal p-p loop. A simple calculation leads us to two important observations: 

1. The second term (equal to A -1 log ~ A -1 ) is negligible as compared to the first one (~ 
A -1 log A). This result can be generalized as follows. If a diagram is reducible in one channel, 
the dominant contributions to its A-derivative will be obtained by deriving only the propagators con- 
necting the irreducible blocks and not the irreducible blocks themselves. This is exactly what we used 
in the derivation of the one-loop flow equations in Section^ 

2. We can replace -Fa (A) in Eq. (1871 by -Fa(O) without changing the result to leading order. In fact, 
integrating the simplified RG equation — d\P(A) — j^F\(0), we obtain indeed the correct result for 
-P(A). This result can be generalized as follows. In the RG equation, one can replace any vertex 
T(ki, ki 1 kz) by its value at £i = £2 = £3 = 0, i.e. the momenta hi can be projected onto the Fermi 
surface. In this sense the dependence of the vertex on the band energies £j is irrelevant. The same is 
true for its dependence on the frequencies. The replacement Fa (A) — > Fa(0) is only justified in the 
differential RG equation Eq. d87l >. but not in the integral expression Eq. d86i . since J d^/^ Fa (0) = 
\og 2 (W/A) ^ -P(A). It is important to note that although the dependence of the vertex on the band 
energies and frequencies is irrelevant in the sense explained above, this dependence does not need to 
be weak at all! In our example, F\(£) depends clearly on £. 



C Linear response 

In this appendix, we investigate the linear response of the system to a weak external perturbation. 

C.l Generalized susceptibilities for superconductivity 

Let the system be coupled to a time dependent, real valued field A(r, t), via an additional term in the 
Hamiltonian 

#ex M = -^A(r,i) (a a (t) +Aji(r)) , (88) 

r 

where 

A A (r)=^/ A (r-r')cv iCrT , (89) 

is the locally defined superconducting order parameter. The index A stands for the internal wave function 
f A (r — r') of the Cooper pair. This term can also be written in Fourier space as 

H mtt t = -V^2(X(-k,t)A A (k)+ h. c. ), (90) 

k 

where 



A A (k) = ^ ikr A A (r) = MP) CpiQc-pt (91) 

r p 
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and the functions A(k) and /a(p) are the Fourier transforms of A(r) and /a(t)- 

The coupling to the field induces in general a non-zero value of (A B (r)) t , where the internal wave 
function fg can be different from /a. The linear response is given by the retarded response function or 
susceptibility XretBA( r ~ v ' ^ ~ *')> v * a 

(A B (r)) t = Jdt'J2xf c ? S BA(r-r',t-t')X(r',t') (92) 

r' 

or, after Fourier transformation in space and time, 

(A B (k)) u = x^Ix(k,w) A(k, w ). (93) 

The global superconducting order parameter is given by A(k = 0). A more convenient quantity is the 
Matsubara response function given by 

X B b C a S (k, v) = dr ( A B (k, t)A a (k) ) , (94) 

where the Bose-Matsubara frequency v is restricted to even multiples of 2n/ ' (3 and the imaginary time 
dependence of operators is given by Eq. M6\ . The retarded response function is obtained by analytic 
continuation 

XW^X^H, (95) 

where 5 is an infmitesimally positive (dissipative) term. 
Equation \9A\ together with Eq. ( I78> . give 



X B B C A S {k) = i rdre-^/^p)/^?') (^(rl^p^rjctp^^ 
p,p' 



(96) 



-rr^y E /b(p)/a(p') G iT (p, k-p,k-p'), (97) 

iE^) [/3^ p ,p' + ^f(p)rf cs (p,p')] iW)7^0, (98) 

p.p' 



where fc = (z/, k), T^ cs (p,p r ) = T(p, k — p, k — p') and /a(p') is the complex conjugate of /a(pO- 

It is argued at the end of Section|3]that the low energy behavior in the weak coupling limit is determined 
by the degrees of freedom close to the Fermi surface. Therefore the dependence of the form factors on 
£ is irrelevant. Let us perform the change of variables p — > £ p , 9 P , where £ is the band energy and 9 
is a suitable (angular) variable. The most elementary superconducting susceptibility is a function of two 
angular variables, which is obtained by putting /s(p) = S(9 P — 9) and /a(p') = S(9 P ' — 9') in Eq. (1981 . 

x BCS (k,9,9') = Jd£ yV 4W.«0iE ^(p (s , e) ,A-p (fl , c) ,A-p (fl ,, eo ) (99) 

PO ,Po 

where J(£, 0) is the Jacobian, such that l/VJ2 p ^ J d0 J d£J(£, 9). 



1 1 We use the conventions F(r) = i J] k e lkr _F(r) and -F(k) = J] r e lkr F(r) for any quantity -F(r). The only exceptions 
are the creation, and annihilation operators. There, we use a pre-factor V~ 1//2 for both transformations. Note that -F^(k) 
denotes the Hermitian conjugate of -F(k) rather than the usual Fourier transform of F^(r). The Fourier transform in (real) 
time is defined as F(t) = /ff e~ iult F(ui) and F(u>) = fdt e iut F(t). 
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Fig. 15 The pattern of charge (spin) currents along the bonds of the square lattice in a charge (spin) flux-phase. 



This susceptibility allows to analyze the superconducting instability without any prejudice on the form 
factor /(p). The natural form factor is obtained by writing x(&, &') m diagonalized form 



x(o,e') = J2xnfn(o)fn(e'), 

n 

where \n an ^ fn are » respectively, eigenvalues and eigenfunctions of the operator 
f{6) -> /d0' x (MW)- 



(100) 



(101) 



If Xn diverges at an energy scale A c , this indicates a (spin-, charge- or superconducting) instability with 
the form factor /„. 

To simplify the calculations, one often selects a certain form factor f(p) and calculates the susceptibility 
for/ A (p) = / s (p)=/(p). 

C.2 Density waves and flux phases 

The effect of a term creating electron-hole pairs is evaluated in the same way as that of a two-particle 
source term. We consider operators of the form 



(102) 



or in Fourier space 

A (q)^E s -/(p) c p 



p(T"p+q(T- 



(103) 



p. a 



The most simple cases are the magnetization (s CT = a) and the charge (s CT = 1), where /(p) = 1, but 
many possibilities with non-trivial p-dependencies exist. Two other examples (for D — 2): if /(p) is a 
function with d x i_ y i -symmetry, q = (tt,tt) and s a = 1 (cr), then the nearest-neighbor-terms in (A(q)) 
yield circular charge (spin) currents flowing around the plaquettes of the square lattice with alternating 
directions (see Fig. I15> . These charge and spin instabilities have been discussed a long time ago in the 
context of the exci tonic insulator |76 1. 

The phase with circulating charge currents is actually identical to the flux-phase 1 77 78 1 and sometimes 
referred to as d-density wave, charge-current wave or orbital antiferromagnetism. Recently this order 
parameter was proposed to compete with d-wave superconductivity and to be responsible for the pseudo 
gap phase of the cuprates |79|. 
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The phase with circulating spin currents is called the spin flux-phase. Other names encountered in 
the literature are "spin current wave" or "spin nematic state" (because it is a state with broken rotational 
symmetry and unbroken time reversal symmetry). The low-temperature thermodynamics of both charge 
and spin flux-phases have been investigated in mean field approximation in Refs. 1 80l l8 II . 

We again considered two different form factors Ja and fs- This allows either to obtain an "angle- 
resolved" susceptibility x s ^ c {q> $j m analogy to Eq. (1991 . or to select the form factor / in advance and 
to choose /a = fs = f ■ The generalized charge- and spin susceptibilities are given by 

xt A (l) = ^ /3 dre-(A B (q,r)(At l (q)+A A (-q))) (104) 

= i J dr fB(p) (T^P 7 ) + f A (p' + q)) • 

p,p' 

■ ^S CT S (T - (cp<T( T ) C P+q<T( T ) C P'+q CT ' C p'<T') ' ( 105 ) 

where q = (y, q). The most prominent case is q = (tt, it, . . .), where A^/^q) can be made Hermitian 
if we choose /(p + q) = /(p). With this choice the spin- and charge susceptibilities can be written in 
perfect analogy with Eq. d98l 



p.p' 



= Tfky £ /s(p) [-f 3 ^' + if W)] <V) Mp')> d07) 

where T s q (p,p') = -XT(jp',p + q,p) and T c q {p,p') = (2 - X)T(p',p + q,p). 

C.3 Flow equations for the susceptibilities 

The RG formalism can be used to calculate the susceptibilities to leading order in the infrared cutoff A. To 
see this, we rewrite Eqs. d98l and (I107i as 

Xl£i(*) = ^E Z S(P)^(P)A(P) (108) 
p 

x^M = -^E<a,W<W^(p)' ( 109 ) 

p 

where the effective field vertices Z are defined as 

zE° A s k ( P ) = Mp) + ^Y,Mp')Dl p k ( P ')rf c k s (p\ P ) (no) 

p> 

< C A » = Mp) ^T.^{v')Dl h q {p')vt q (p'M- (in) 

p' 

The vertex function Y BCS in Eq. (II 10> can be expressed by an infinite series in the p-p irreducible part 
I pp . Formally 

ZE° S = f B + f B D pp T BCS 

= Jb + f B D pp I pp + f B D pp I pp D pp I pp + ■■■ 

= f B + Z§ cs D pp I pp (112) 
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To get the flow equation, we take the derivative with respect to A of Eq. (II 12t . It gives three terms 

Z BCS = z BCS D pp I pp + Z BCS D PP I PP + Z BCS D PP I PP . (113) 

The last term of Eq. J 1 1 3i is neglected for the same reason as in Eqs. (1421 . The remaining equation is 
iterated to give 

gBCS _ 7BCS Jjppjpp 1 rrBCS jjppjpp TjPPTPP I 7BCS TjPP WP TjPP TPP TjPP TPP I . . . 

B B B B 

= Z BCS D PP T BCS (114) 
This is the RG equation for the field-vertex Z§ cs . The form factor Jb enters as an initial condition 

zfs?, fc (pw*(p). 

It is now easy to obtain a RG equation for x BCS as well. From Eqs. dl08> and il 14i 

■BCS ?BCS npp f , 7 BCS riPP f 

X = Z B D JA + ^B L ' JA 

= Z BCS D PP T BCS D pp f A + Z^ cs D pp f A 

= zE cs D pp (f A + r BCS D pp f A ) 

= Z BCS D PP Z BCS (115) 

In the last equation, we have introduced the function Z B< ]f k (p), which is the same as Z§ c A s k (p), except 
that the initial condition /s(p) is changed into f A {p). 

Written out, the RG equations for the pairing susceptibility are 

= ^zZ Z B%>)D^(p)Z%Z(P) (HQ 
p 

Z^Up) = ^J2z^ S k (p')Df tk (p')T BC k s (p\ P )^ (117) 
p' 

where o = A or B, and T BC j[ s (p' ,p) — T A (p' , k — p' , k — p) . Once the vertex function Y BCS is known from 
the one-loop RG equation, one can solve the linear equation M 17i for the initial condition Z B ? s k (p) = 
/o(p) and finally one can integrate Eq. il 16ft to obtain the susceptibility. 

All the steps can be repeated for the charge- and spin susceptibilities, just by replacing the quantities 
Z BCS , T BCS , I pp by Z s / C , T s / C , l s ' c and the p-p propagator D pp by -D ph . Thus, 

p 

K'x,m = 4K:>KM r S(p'^ (11 9 ) 

where o = A or B, T% q {p',p) = -XT A (p',p+ q,p) and T\ q {p',p) = (2 - X)T A (p',p + q,p). 

D Functional integral formulation 

In the functional integral formulation the annihilation- and creation operators c kcr and are replaced by 
anti-commuting Grassmann fields VVr.k.o-) an d ^(i-.k.o-)- They depend on the "imaginary time" variable 
r 6 [0, 0] in addition to the labels of the single-particle states, and they satisfy anti-periodic boundary 
conditions in the r variable. 

The Hamilton operator is transformed into the action 

S[lp] =~ J\t fX)^(r,k,«r)(ar-M)V'(r,k, < r)+fl"[T,V]^ » (120) 
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where H[t, ip] is obtained from the Hamiltonian by the replacements — > V(r,k,<r) ar, d — > ip( Tj u. !(T ). 
It is useful to take the Fourier transform 

4><rk=0- 1/2 / dre^W.k..), (121) 
Jo 

where k = (ko, k) contains the Matsubara frequency fco- It allows to write the action as 

S[V] = 5oM-W / M, (122) 

with 

s o[ip} = X! ^fc( ifc o - &)ip<rk, (123) 
where £k = ek — /i and 

W W\ = \~py 9{^U^2,^3)^24>k 1 ^Jk 2 a'^k 3 a'^k 1 +k 2 -k 3 a- (124) 

The partition function is given by the functional integral 

Z = [ V^e sw , (125) 



where Dip is a short-hand notation for J\ kcr dipkadipka- The one- and two-particle Green's functions are 
given by 

G(fc) = -(V^fer) (126) 



G" TCr (fcl,fc2,fc 3 ) = /3V (■>pk 1 a'tpk2<T''4>k 3 a>'<Pk 1 +k2-k 3 cr)- (127) 

where the average of a Grassmann expression 0[^] is defined as 

W]} = i / DVe^O^]- (128) 



It is convenient to introduce the partition function with source term 

Z[rj\ = J dnc[*P] e - w ^ + ^ + ^' r >\ (129) 

where we used the short-hand notation (x, tp) := X) CT fc Xaki>ak and the normalized Gaussian measure is 
defined by 

^ : = JVteVW (130) 
The connected part of the average of any Grassmann monomial is obtained as a functional derivative 1 60 1 

5 2n log Z [17] 



(V>1 • • -Ipnlpn+t ■ ■■1p2n)c = 



where we have written ipi instead of tpkio-i ■ 



(131) 

r)=0 
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